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1 Introduction 



The study of event shape distributions in e + e~ annihilation has proved to be a valuable 
testing-ground for QCD. For an infrared-safe shape variable, i.e. one that is insensitive 
to the splitting of a final-state momentum into collinear momenta and to the emission 
of soft momenta, the distribution can be calculated order-by-order in perturbation theory. 
However, one finds that higher-order contributions are enhanced by large logarithmic factors 
near the two-jet region, where most of the data reside. Hence for a detailed quantitative 
understanding one should resum as many of these terms as possible to all orders. This 
has been achieved for a limited number of shape variables, namely those shown to satisfy 
an exponentiation property, to be specified in more detail later. These variables include 
the thrust [0, heavy jet mass 0, jet broadening [[|, energy-energy correlation ||] and 
differential two-jet rate f5|. 

In the present paper we show that the C-parameter || shape variable also has the 
required property of exponentiation, so that we are able to carry out the resummation of 
large logarithms to the same accuracy as for those mentioned above. The exponentiation 
of soft gluon contributions to QCD matrix elements is a general phenomenon; the critical 
property for event shapes is that the corresponding phase space should factorize with 
sufficient accuracy to maintain exponentiation of at least the leading and next-to-leading 
logarithms. In the case of the C-parameter, we are able to show this by exploiting a simple 
connection between the C-parameter and the thrust in the dominant phase-space regions, 
namely those where C is small. 

In Sect. |2| we review the definition and kinematics of the C-parameter, and the rela- 
tionship between it and the thrust in the small-C region. Sect. |3] presents the fixed-order 
predictions for the C-parameter distribution up to order a|, together with the large loga- 
rithmic terms that appear at small C. The resummation of these terms, to next-to-leading 
logarithmic accuracy, is performed in Sect. |4] using the connection with the thrust. 

In order to describe the C-parameter distribution over the widest possible range, one 
should match the resummed results to those at fixed order outside the two-jet region. A 
suitable matching procedure, chosen from those proposed in Ref. is outlined in Sect. ^. 

A further improvement in the prediction can be made by estimating non-perturbative 
effects, which have been found to be substantial for many shape variables at present en- 
ergies. This is discussed in Sect. |6], where we argue that the effects on the C-parameter 
distribution can again be related to those for the thrust, via their connection in the soft 
region. 

Finally in Sect. [7| we make a preliminary comparison with experimental data and discuss 
the results. 
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2 Kinematics of the C-parameter 



The C-parameter was initially defined as || 

C = 3(AiA 2 + A 2 A 3 + A 3 Ai) (1) 
where \ a (0 < A Q < 1, J2a ^« = 1) are t ne eigenvalues of the linearized momentum tensor 



_ Pi Pi /\Pi\ ( . , 



the sums being over all final-state particles. Neglecting particle masses, we may express C 
in terms of invariants as 

where Q M is the total four-momentum 

Q" = £rf- ( 4 ) 

Introducing the cm. energy fractions X{ = 2{pi ■ Q)/Q 2 , we can write Eq. (g) as 

C — I sin 2 0^- . (5) 

The kinematic range is < C < 1, with C = for a perfectly two-jet-like final state and 
1 for an isotropic and acoplanar distribution of final-state momenta. In fact the maximal 
value C = 1 can only be achieved when there are four or more final-state particles. Planar 
events fill up the kinematic region C < |. For three particles, the maximum value C — | 
corresponds to the symmetric configuration x% — x% — x% — |. 

2.1 C-parameter and thrust at small C 

Let us consider the definition (ffl) of the C-parameter and assume the ordering of eigenvalues 
Ai > A 2 > A3. It follows that in the small-C limit the maximum eigenvalue approaches 
unity, Ai — > 1, whilst A 2 and A 3 are parametrically small, A 2 ~ A 3 ~ 0(1 — Ai). Using the 
normalization condition Y> a \ a = 1, we thus obtain 

C = 3{(1 - AO - i[(l - AO 2 + A 2 + A 2 ]} ~ 3(1 - A0[1 + 0(1 - AO] , (6) 

or, equivalent ly, 

C = 3(1 -A x )[l + 0(C)] , (C-0). (7) 
We can now relate the eigenvalue Ai to the thrust. The definition of the thrust T is [[F] 

T = Max S,t|Prn| , (8) 



where the maximum is with respect to the unit 3-vector n. At the maximum n = 
the thrust axis. Denoting by ri\ the eigenvector of the linearized momentum tensor (|2|) 
corresponding to the maximum eigenvalue Ai, by definition we have 

1 „ \p, ■ nA 2 E,-|p,--ni| m . , 

Ai = — r ^» ■ r ^ ' , < r , (9) 



^j\Pj\ \p, 

where the last inequality comes from Eq. (§). From Eq. (0) we thus obtain 

C > 3(1 -T)[l + 0(C)] , (C->0). (10) 

Equation ( ]TD| ) implies that in the small-C limit the thrust approaches unity. The oppo- 
site is also true. In order to show this, let us consider the separation of the final state into 
two hemispheres S+ and S- by the plane orthogonal to the thrust axis tit- We can write 
Eq. (D as 

C = - (>:', ; .r,.r, sin 2 % + Tf^XiXj sin 2 6^) , (11) 

where Y^'ij (Sfj) denotes the sum over all final-state pairs belonging to the same (opposite) 
hemisphere (s). Then, applying the identity sin 2 9y = 2(1 — cos 9y) — (1 — cos 6^) 2 (sin 2 6^ = 
2(1 + cos 6ij) — (1 + cos%) 2 ) whenever the particles i and j belong to the same (opposite) 
hemisphere (s) and using energy- momentum conservation, it is straightforward to recast 
Eq. ( |TT1) in the form 

C = 3 [ 2(w + + w-) - (w+ - W-) 2 ] - A T , (12) 
where w + and W- are the rescaled invariant masses-squared of the two hemispheres, 

^HsA ' <i3) 

and A T is defined by 

A T = - (£■ jXiXjil - cos%) 2 + £".2^(1 + cos%) 2 ) . (14) 



Note that A^ is positive definite and that the expression in the square bracket in Eq. (|Tj) 
is exactly equal to 1 — T 2 |IJ. Therefore from Eq. (|T^p we obtain 

C = 3(1 - T 2 ) - A T < 3(1 - T 2 ) < 6(1 - T) . (15) 



Thus, combining the inequalities ( |TUD and (IB), we see that for sufficiently small C (i.e., 
neglecting corrections of order C 2 ) we have 

3(1 — T) < C < 6(1 — T) . (16) 
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3 The C-parameter in perturbation theory 



The C-parameter is manifestly insensitive to the splitting of a final-state momentum into 
collinear momenta (9ij = 0) and to the emission of soft momenta (xi = 0). It follows that 
the C-parameter distribution can be computed in QCD perturbation theory as a power 
series expansion in the strong coupling as. For (7^0, the predicted distribution has the 
general form 

~ = a s A(C) + a 2 s B(C) + 0(a s ) , (17) 

(To 

where a$ = as/27r and we normalize to the Born cross section <To, as was done in Refs. MM- 
The first-order distribution A(C) is shown by the dashed curve in Fig. |l|. Notice that it 
diverges as C — > 0; in fact at small C one finds 



A(C) = 4C F i 



+ C(lnC) . 



At larger values of C, A(C) smoothly approaches a finite value at the three-parton upper 
limit C = §: 

Ml) = 2if ^ = 7 - 6433 • ( 19 ) 

It will be useful to define the fraction of events with C-parameter values less than C: 

R (C) = f° — ^ = 1 + a s Ri{C) + a 2 s R 2 {C) + 0(a|) (20) 
Jo a t aC 

where a t is the total cross section. Then since -R(C max ) = 1 and <7 t /<J = l + ^Cpas + O^ag), 
we have 



RAC) = - f 4 A(C)dC 
Jc 

i - 

R 2 (C) = - B(C)dC +|C F ( 4 A(C)dC . (21) 

J C J c 



Introducing 



we find as C — > 



In ( | ) (22) 



i^i(C) = -2C F (L 2 - \L - In 2 + |) + O(ClnC) . (23) 

Notice that in first order we obtain up to two large logarithmic factors L at small C, 
corresponding to the emission of one collinear and/or soft gluon. Note also that by rescaling 
C by a factor of 6 in the definition (p2] ) we absorb all logarithms in Eq. ( |23| ) into L. In Sect. f| 
we shall see that this is a consequence of the relationship fll5|) between the C-parameter 
and the thrust at small C. 

Upon adding the second-order contribution B(C), we obtain the solid curve in Fig. [TJ, 
which exhibits a number of new features IBfl: the behaviour at C —>■ is modified, in fact 
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Figure 1: Fixed-order predictions of the C-parameter distribution for as = 0.12. 
Dashed: 0{a s ). Solid: 0(a|). 

becoming strongly negative divergent; the distribution remains finite as C — > | but there 
is a new divergence as C — > f + ; and finally the distribution vanishes smoothly as C — > 1. 

The new features at large C arise from the abrupt change in the 0(a,s) contribution 
A(C) at C — |, see Eq. ([19]). They have been discussed in Ref. |T0] and the resummation of 



large terms of the type ln(C — |) in this region will be considered in a separate paper [[11 
In the present paper we concentrate on terms that diverge as C — > 0, and the matching 
of the resummed and fixed-order predictions in the region < C < f , where the latter 
remains smooth. 

In the region of small C, the emission of two collinear and/or soft gluons can yield up 
to four powers of L = ln(6/C). Thus, writing the logarithmic dependence explicitly, the 
second-order contribution R 2 to the event fraction R(C) defined in Eq. (|20|) has the form 

R 2 {C)= R2mL m + D 2 {C) , (24) 

where D 2 — >• as C — >• 0. We have computed the coefficients R 2 4, R 2 3, R 22 analytically, and 
find that their values, together with Eq. ([23), are consistent with the following exponenti- 
ating form for R(C): 

(co \ / oo n+1 \ oo 

1 + £ C n a n s exp £ Gn m a n s L m + ^ a n s D n (C) . (25) 
n=l / \ra=l m=l / n=l 

Here C n and G nm are C-independent coefficients and the remainder functions D n (C) vanish 
as C -> 0. 



In the next section we argue that the exponentiation formula ( p5|) is actually valid to 
all orders in perturbation theory. Therefore it is convenient to express the coefficients in 
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G\2 

G22 
G23 



+C F (4n 2 - 15)/6 

+3C F 

-2C F 

-C f [A8tt 2 C f + (169 - \2tx 2 )C a - 22iV/]/36 
-C F (llC A -2N f )/3 



Table 1: Coefficients C\ and G nm in Eq. ([251). 



Eq. ( p4|) in terms of those in Eq. (^), as follows 



-R20 
-R23 



C2 , -R21 — C21 + C\G\\ , 
G23 + C12G11 , -R24 



22 



G22 4- i^ii Ci^* 



12 , 



lp2 
2^12 



The coefficients which determine i?2m for m > 2 are as given in Table 1. 



(26) 



The coefficients C2 and G21 which enter into R20 and R21 are not known analytically but 
can be fitted to numerical data from the Monte Carlo matrix element evaluation program 
EVENT @| . From Eqs. (|T]) and (||) we have, for N f = 5, 



- f dCB{C) = R 2 {C) + lC F R 1 {C) 



(27) 



c -» 



3.556L 4 - 20.889L 3 - 36.778L 2 + (G 21 + 29.758)L + C 2 + 10.879 . 



By fitting we find 

C 2 = 76.5 ± 2.9 , G21 = 63.4 ± 6.0 . (28) 
Then Eq. (^) gives a good description of the matrix element, as illustrated in Fig. |2|. We 



checked that similar results are obtained with the program EVENT2 [12 



4 Resummation for C 







As explained in the Introduction, in order to carry out the resummation of the logC 
contributions for C — > 0, we exploit the connection between the C-parameter and the 
thrust in the two-jet region, discussed in Sec. 0. However, to derive the relevant relation 
between C and 1 — T, we cannot simply rely on the kinematics but rather we have to 
consider QCD dynamics. 

As discussed in detail in Ref. [fjj, the structure of the soft and collinear singularities in 
the multiparton QCD matrix elements can be described in terms of a branching process. In 
the large-T region and to next-to-leading logarithmic accuracy, this process takes place in an 
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C 

Figure 2: Second-order prediction of the C-parameter distribution for C — > 0. 
Points: EVENT Monte Carlo. Curve: Eqs. (£7}J(8fy. 

angular-ordered region around the thrust axis. Because of this angular ordering (enforced 
by QCD dynamics), we can introduce in Eq. (|Tl|) the approximation sin 2 Oij ~ 2(1 — cos#y) 
(sin 2 ^- ~ 2(1 + cosOij)) in the terms appearing in Yhj (J2'lj)- Equivalently, we can say 
that when T — > 1 the phase-space region contributing to Ay in Eq. (|15|) is sub dominant 
with respect to that contributing to 3(1 — T 2 ) ~ 6(1 —T). It follows that to next-to-leading 
logarithmic accuracy at small C we have 

C~6(l-T). (29) 

In Ref. [0 the perturbative expression for Rt{t), the fraction of events with thrust 
larger than T = 1 — r, was written in the following form: 

R t (t = 1 - T) = K T (a s )Z T (T, a s ) + D T {r, a s ) . (30) 

By definition the functions Kt(ocs), Dt(t,cxs) and Et(t, as) are power series expansions 
in as whose coefficients are respectively constant in r, vanishing for r — > and polynomials 
in L = In 1/r. The perturbative contributions to Rt{t) which are logarithmically enhanced 
in the two-jet region r = 1 — T — > are thus embedded in Ey. At small r it becomes 
important to resum the series of large logarithms in Sy. 

By naive counting of logarithms Rt, and hence Et, contain terms of the type a^L™ 
with m < 2n. However, by explicitly performing the resummation of the leading logarithms 
(those with n < m < 2n), it was shown [Q] that they exponentiate. The word exponentiation 
means that the terms a s L m with m > n + 1 are absent from lnR? (lnE^), although they 
do appear in R T itself. 

The result in Eq. (p9| ) allows us to obtain the leading and next-to-leading logC con- 
tributions to the C-parameter distribution simply by replacing ln(l/r) by L = ln(6/C) in 
the corresponding formula for the thrust distribution. It follows that the function R(C) 
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defined in Eq. ( pO|) has the same form as Eq. fl30|) , that is, 

R(C) = K(a s )E(C, a s ) + D(C, a s ) , (31) 

with 

oo n+1 

ln£(e> s ) = J2J2 G nma n s L m (32) 

n=l m=l 

= Lgi(a s L) + g 2 (a s L) + a s g^{a s L) H . 

The function ft resums all the leading logarithmic contributions a"gL n+l , while g 2 contains 
the next-to-leading logarithmic terms atgL n , and ft etc. represent the remaining subdom- 
inant corrections otgL m with < m < n. All the functions ft vanish at L = since they 
resum terms with m > 0. 

Because of the result in Eq. (f29|), the functions ft and ft in Eq. (|32"[) are the same as 
those for the thrust distribution, given in Ref. provided we identify L with m(6/C). 
This relation between lnE(C) and lnS^(r) breaks down at subdominant orders, i.e. the 
functions ft for i > 2 are different in the two cases. Similarly the non- logarithmic coefficient 
functions K and and the remainder functions D and are different, as may be seen 
by comparing the fixed-order results in Sec. |3| with the corresponding ones in Ref. [|IJ . 



5 Matching with fixed order 

To combine the fixed-order and resummed predictions of the previous two sections without 
any double counting, we may adopt any of the matching procedures proposed in Ref. 
In the so-called log-i? matching scheme, for example, one writes the next-to-leading order 
expression fl2"U|) in the equivalent form 

\nR(C) = a s Ri{C) + a 2 s {R 2 {C) - ^(C)] 2 } + 0(af) (33) 

and then replaces the computed logarithmic terms by the corresponding resummed expres- 
sions, to obtain 

\nR(C) = L gi (asL)+g 2 (asL)+as{Ri(C)-G u L-G 12 L 2 } 

+a 2 s {R 2 (C) - \{Ri{C)} 2 - G 22 L 2 - G 23 L 3 } . (34) 

An advantage of this scheme is that, since all terms are exponentiated, it is not necessary 
to separate out the coefficient and remainder functions K and D in Eq. ( |31~D explicitly. A 
potential disadvantage is that the resummed terms may dominate over the fixed-order ones 
far from the two-jet region, if the latter become very small there. This should not be a 
severe problem for the C-parameter distribution, since the first-order prediction does not 
vanish at the three-parton boundary, Eq. ([19]). Thus we expect the log-i? prescription ([34]) 
to be satisfactory at high energies, as long as matching is limited to the region C < |. 
Note, however, that the sensitivity to the matching procedure can increase at low energies, 
and therefore other schemes should also be tried. 
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6 Non-perturbative power corrections 



Even at high energies (Q^,Mz), it has been found that substantial non-perturbative 
'hadronization corrections' need to be applied to perturbative predictions concerning event 



shape variables. Traditionally these have been estimated [13,1^] using Monte Carlo hadroniza- 
tion models [|I^,|16[]. Recently, an alternative 'renormalon' or 'dispersive' method for esti- 
mating them has been developed |I7| - |25| . According to this approach, the dominant non- 



perturbative correction to the thrust distribution, at thrust values that are not too close to 
the kinematic boundaries, should correspond to a simple shift in the resummed perturbative 
prediction, of the form |24 

R T (r) = Rr\r-4A 1 /Q) , (35) 

where i?jf rt is the function called Rt in Eq. fl3~0|) and A\ is a non-perturbative parameter 
to be determined experimentally 

Because of the close connection (|16D between the C-parameter and thrust, the same 
type of result holds for the C-parameter distribution. Furthermore, the assumption that 
the leading non-perturbative effect is associated with a universal effective strong coupling at 



low scales 0] leads to a relationship between the shifts in the two distributions. Consider 
the emission of a single soft gluon at angle 9 with energy fraction x. We have from Eq. (|^) a 
contribution to C of |x sin 2 9, whereas the contribution to r is \x min{l — cos 9,1 + cos 9}. 
As expected, the ratio is between 3 (at 9 = tt/2) and 6 (in the collinear regions, 9 = 0,n). 
Integrating at a fixed small value of the gluon transverse momentum k t , we obtain 



SC \ 2Cp ^ ri dx r J9_ _ ^ p , [ f * sin 2 9 



St n Jo x Jo sm9 UKt I* min{l - cos0, 1 + cos0} 



giving ]T7| 



(36) 
(37) 



6C_ _ 3tt 
57 ~ T ' 

This is smaller than the factor of 6 in Eq. (|29|), because the non-perturbative effect comes 
from the whole soft region of gluon emission, whereas the perturbative logarithmic en- 
hancement comes from the collinear region. As in the case of the thrust distribution, the 
soft-gluon contribution exponentiates and, taking into account Eq. (|37D, we have 

R(C) = R pcvt {C - QtxA 1 /Q) . (38) 



In Ref. h was pointed out that higher-loop effects can alter the coefficients of the 
power corrections to event shapes predicted by the dispersive approach. However, recent 
studies of two-loop contributions |26j have shown that Eqs. ( j35|) and ( |3~8D remain valid, 
provided that the relationship between the non-perturbative constant A\ and the effective 
coupling at low scales is renormalized by a common overall 'Milan factor' M. ~ 1.8. 



7 Comparison with data 

At present the experimental data on the C-parameter distribution are limited to high 
energies, Q = Mz [14,27] and above [28]. Figure |3] shows the LEP1 data [^j]] together with 
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various theoretical calculations, all assuming as(Mz) = 0.118. The solid curve represents 
Eq. (^) with the best-fit value of A\ = 0.24 GeV, which gives very good agreement 
throughout the region 0.09 < C < 0.75. For the function R pert the resummed expression 
with log-i? matching, Eq. fl34|), was used. The dashed curve shows the effect of neglecting 
the non-perturbative shift Ax, while the dot-dashed curve shows the fixed-order prediction 
(|T7j), with neither resummation nor non-perturbative correction. 



The corresponding value of A\ obtained from fitting the thrust distribution EM is A\ 



0.22 GeV. Thus the results obtained from the C-parameter and thrust data are consistent 
within 10%. Taking into account the Milan enhancement factor M. in the coefficient of 
the power correction p6| , these results correspond to a value olq ~ 0.35 for the non- 
perturbative parameter d (the mean value of the effective strong coupling at scales below 
2 GeV) introduced in Refs. MM. 



The fitted values of as and A\ are not well constrained separately by the existing C- 
parameter data, owing to the lack of measurements of this quantity at lower energies. A 
comprehensive analysis over a wide range of energies, similar to those presented in Refs. 
T| for other shape variables, would therefore be very useful. 
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